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ABSTRACT: The spectrum of anomalous dimensions of twist st(2) operators in M = 4 
SYM has an intriguing feature in low twist 2 or 3. The anomalous dimension of the lowest 
state, dual a folded string on AdSs x S 5 , can be computed by Bethe Ansatz at 3, 4 loops 
respectively as a simple closed function of the Lorentz spin. This feature is apparently 
lost at higher twist. We propose sum rules for the excited anomalous dimensions where 
closed expressions can still be provided, even at higher twist. We present several explicit 
three loop examples. Many structural regularities can be observed leading to closed ex- 
pressions which depend parametrically both on the spin and the twist. They allow to 
compute the subleading term in the logarithmic large spin expansion of the sum rules as 
a compact simple function of the twist, in analogy with the recent results by Freyhult, Rej 
and Staudacher in arXiv:0712.2743 [hep-th]. 
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1. Introduction 

The long range Bethe Ansatz solution of the mixing problem in M = 4 SYM allows to 
compute multi-loop anomalous dimensions of various composite operators in a very effi- 
cient way [p. The main and unique obstacle is the celebrated wrapping problem [|p setting 
an upper bound on the achievable order in the loop expansion. Within this bound, the 
exact perturbative anomalous dimensions are recovered without approximation. 

In some cases, we areinterested in parametric classes of operators where we would 
like to compute the spectrum as a closed function of the characterizing parameters (Lorentz 
spin, operator length, etc.). This is a much more difficult question than merely asking the 
value of the anomalous dimensions at a specific point in the parameter space. 

A positive answer beyond one-loop is not known in the general case, but is available 
for certain specific classes of operators where additional insight saves the day. A very 
simple example is that of the generalized gluon condensate 

Ol = Tr T L + higher order mixing, (1.1) 

where T is a component of the self -dual field strength [^|, ||, ^, |p. In this case, it is possible 
to derive compact five loop expressions for the anomalous dimension as a closed function 
of the parameter L [0J. 

A much wider and more interesting class is that of quasi-conformal operators [0]. 
They span an integrable sector of QCD and its various supersymmetric generalizations, 
including of course TV = 4 SYM. In the cases relevant to our discussion, they take the form 

O n = J2 c nu ..., nL Tr{D^X---D nL X), (1.2) 

niH h"L 

where D is a light-cone projected covariant derivative and X can be an elementary scalar 
(ip) or a so-called good component of the gaugino (A) or gauge (A) fields. In this context, 
L is the twist of the operator and N is the total Lorentz spin. 

The case when X is a scalar identifies the 5 [(2) sector which is closed at all orders in the 
gauge coupling [B]. The gaugino case describes special operators in the purely fermionic 
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closed s((2|l) subsector which evolve autonomously under dilatations. Finally, the gauge 
operators close at one-loop and have a complicated pattern of higher order mixing. 

In general, the one-loop classification of states is made simpler by the underlying 
collinear s[(2) algebra under which the charges of <p, A, A are s = |, 1, 1 respectively 
The twist L states belong to [s]® L and can be decomposed in irreducible infinite dimen- 
sional sl(2) modules. The twist-2 case is very special since supersymmetry links together 
the three physical values of s [ |10| , [lip. One has a single supermultiplets and a universal 
anomalous dimension 7 un iv(-^) describing (with trivial shifts in the Lorentz spin) all the 
states. This anomalous dimension can be computed at 3 loops as a closed function of N 
by invoking the Kotikov, Lipatov, Onishchenko and Velizhanin (KLOV) principle leading 
to a simple Ansatz in terms of nested harmonic sums [12]. 

The case of twist-3 is more complicated. There are in principle three distinct series 
of modules associated to each elementary field [jl3|]. Detailed results, including a nice 
application of superconformal symmetry, are described in [14, 15, [TJ] for the scalar sector, 
[ p7| ] for the gaugino sector, and [jl8|, [l^] for the gluon sector. 

A curious fact is that the ground state (i.e. that with smallest anomalous dimension) 
admits again closed expressions in all sectors based on suitably generalized KLOV-like 
principles. As soon as one moves to higher twist, even in the simplest case of X = (p, 
it is easy to check that no simple closed expressions describe the ground state whose 
anomalous dimension is irrational. Therefore, the following basic questions naturally 
arise: 

1. Why are twist-2 and 3 so special ? 

2. Can we generalize the twist-2 and 3 closed expressions to higher twists ? 

The aim of this paper is precisely that of answering the above questions. We shall 
show that, at higher twist, it is possible to consider sums of (powers of) anomalous di- 
mensions of the ground and excited states. For these combinations, we shall provide 
quite simple closed formulas as well as compact sum rules which are parametric in both 
N and the twist L 

These expressions are a simple hidden constraint on the anomalous dimension whose 
precise nature is not clear. In particular, they suggest similar sum rules at strong coupling 



for the dual string states discussed for instance in [|2(], ^ |22|, g^, [24| ]. 

The plan of the paper is the following. In Sec. @ we recall some basic facts concerning 
the XXX_ S integrable spin chain and provide several explicit examples in Sec. (0). The 
outcome of this analysis is summarized in Sec. @. In Sec. (||) we recall a few important 
properties of the large spin expansion of twist anomalous dimensions. In Sec. @ we 
propose linear sum rules for the singlet anomalous dimensions of various twist operators. 
In Sec. @ we present three loop results for these sum rules in the scalar sector. These 
results suggest a twist-dependent conjecture formulated in Sec. @. Various checks are 
performed in Sec. @. The subleading corrections at large spin are computed in Sec. (|lC|). 
Quadratic sum rules are proposed in Sec. (11), elaborated in Sec. dl2| ) and checked in 
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Sec. ([l3|). A few results for cubic sum rules are finally presented in Sec. (14). Various 
Appendices are devoted to some technical results. 



2. The integrable XXX_ S chain 

In this section, we briefly recall a few basic facts about the integrable XXX- S spin chain. 
It describes the one-loop mixing of twist-L quasipartonic operators built with elementary 
collinear conformal spin s fields. A nice and accessible recent review on this standard 



material is [|2g]. Our presentation is brief and just sets up the language for the later dis- 
cussion. 

2.1 Basic facts 

The XXX_ S chain is a quantum spin chain with s((2, M) symmetry. Each site carries the 
infinite dimensional [s] representation of the 5L(2,M) collinear conformal group. The 
decomposition rule for the tensor products of this representation reads 

oo 

[*]»[*] = ©[* + "], (2-1) 

n=0 

and can be used to analyze the states of a chain with L sites. The local spin chain integrable 
Hamiltonian reads 

L 

,n+l> 

(2.2) 

n=l 

with 

#n,n+l = ^(Jn,n+l) " ^(2 *), 1p{z) = — log T(z). (2.3) 

The quantity J n ,n+i is the spin of the two-site states 

(S n + 5„_|_i) 2 = J n ,n+1 (Jn,n+1 ~ !)• (2-4) 
This Hamiltonian is integrable and can be studied by Bethe Ansatz or Baxter techniques. 

2.2 Relation with M = 4 SYM operators and anomalous dimensions 

The XXX- S Hamiltonian is the mixing matrix for the one-loop renormalization of planar 
composite operators of the form 

<D L {N)= c ni _ nL Tr{D^X---D^X}, (2.5) 

niH \-n L —N 

where 5 = 1/2, 1, 3/2 for the physical cases X = ip, A, A. A straightforward application of 
Eq. Q leads to 



1*]® L = © 9l(N) 



N=0 



l + N 
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which tells that there are highest weight states with Lorentz spin N. These are in 1- 

1 correspondence with the non-trivial (i.e. without roots at infinity) solutions of the Bethe 
Ansatz equations (BAE) 

u k + is \ L _ u k - u 3 - i ^ 27 ^ 
^Uk-isJ l = A Uk-Uj + i' 

A short calculation also provides the energy and momentum charges 

£ = E^T' e ip =f\^^-. (2.8) 
These energies are the one-loop anomalous dimensions of the above operators, say 

7 = ^-^, (2-9) 

where A is the 't Hooft large N c coupling. 

Clearly, it is complicated to enumerate the full set of solutions to the BAE. A useful 
alternative method is the Baxter approach described in the next section. 

2.3 The Baxter equation for the XXX_ S chain 

An alternative approach to the solution of the BAE is based on the Baxter approach [^6|]. 
The main tool is the Baxter operator whose eigenvalues Q(u) obey a relatively simple 
functional equation. If Q(u) is assumed to be a polynomial, then the Baxter equation is 
equivalent to the algebraic Bethe Ansatz equations for its roots to be identified with the 
Bethe roots. A more general discussion can be found in p5| , pffi |. 

In practice, the Baxter approach in the present case is quite simply stated. One intro- 
duces the Baxter function which is the minimal polynomial with roots equal to the Bethe 
roots 

N 

Q(u) = Y\(u - u k ). (2.10) 

fc=i 

The BA equations are equivalent to the Baxter equation that we write for general confor- 
mal spin s (although we shall be mainly interested in the case s = 1/2) 

(u + i s) L Q(u + i) + (u-i s) L Q{u -i) = t(u) Q(u). (2.11) 

where 

t(u) = 2u L + q 2 u L - 2 + q 3 u L ^ + ■ ■ ■ + q L , (2.12) 

q 2 = -(N + Ls){N + Ls - 1) + Ls{s - 1). (2.13) 

The quantities 93, • • • , <7l have the meaning of quantum numbers. They must be obtained 
by consistence of the Baxter equation and the assumption of a polynomial Baxter function. 
Once Q is found, the energy and momentum can be written in terms of Q as 

E = i [(log Q(u))'] + _i s s , e lP = §p^f (2-14) 
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3. The XXX_ S chain at twist L = 2, . . . , 5 

In order to set the stage for the later discussion and introduction of sum rules, we now 
illustrate very explicitly the structure of the Baxter equation at various small twists L = 
2, . . . , 5. Our results will be generically valid positive values s > of the (quantized) 
conformal spin. 

3.1 Twist 2 

For L = 2 and Lorentz spin N, the Baxter equation is 

(it + i sf Q(u + i) + {u-i sf Q(u - i) = t(u) Q(u), (3.1) 

with 

t(u) = 2u 2 + q 2 , (3.2) 

q 2 = N-N 2 -4Ns-2s 2 . (3.3) 

In this case, there are no additional quantum numbers. This follows from the trivial mul- 
tiplicities in 

oo 

[s]®[s]= ®[2s + N]. (3.4) 

N=0 



The Baxter polynomial with degree N (even or odd) is [ ]28| ] 

-N N + 4s - 1 s - iu 



Q{u) = 3 F 2 
The Baxter polynomial has parity 



2s 2s 



1 . (3.5) 



Q(-u) = (-If Q(u) — > e* p = (-l) N . (3.6) 

The energy is 

E = i [(log Q(u))'] +;* = 4 [ip(N + 2 s) - i>(2 s)} , (3.7) 

where 

rf>(z) = £ log r(z). (3.8) 

3.2 Twist 3 

For L = 3 and spin JV, the Baxter equation is 

(u + i sf Q(u + i) + {u-i s) 3 Q(u - i) = t(u) Q{u), (3.9) 

with 

t(u) = 2 u 3 + q 2 u + q 3 , (3.10) 
q 2 = N - N 2 - 6Ns - 6s 2 . (3.11) 
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In this case, there is an additional quantum number. This is related to the multiplicities in 

oo oo 

[a]®[a]®[s]= © [3s + n 1 +n 2 }= @(N + l)[3s + N}. (3.12) 

Indeed, looking for a polynomial solution to the Baxter equation we find the condition 

P(q 3 ) = 0, degP = N+l. (3.13) 
Let us analyze the solutions for real 5 > 0. 
even N 



For even N, the polynomial P reads 

= ft fl(<fe), R{q 5 ) = R{-q 5 ). 



(3.14) 



Hence we can have q% = or q$ = ±q* for some values of q* ^ 0. The Baxter function 
associated with q 3 = is even and is associated with a zero momentum state which turns 
out to be the ground state. 

The explicit form of Q(u) is known and reads 



\ 2 + S 2 + S 2 + S 



Q(-u). 



(3.15) 



To compute the energy, it is convenient to relate Q to the Wilson polynomials 

W n (u 2 ,a,b,c, d) _ I 
(a + b) n (a + c) n (a + d) n 4 3 ' 

We find apart from a trivial scaling 



n n + a + b + c + d— 1 a + iu a — iu 
a+b a+c a + d 



Q(lt) = Wjtya U ,5,5,5, - . 



Using the following formula from the Appendix B of [29] 

iy- W n (u , a, a, c, d) 
= ip{n + a + c) — tp(a + c) + tp(n + a + d) — tp(a + d) 
and the fact that W is invariant under permutations of a, b, c, d, we find the result 



E = i [(log Q(u)) 

'N 



n+is 

i s 



+ 25 - V(2s) + ^ 



j\r 1 



^ 5 + 



(3.16) 
(3.17) 

(3.18) 
(3.19) 

(3.20) 
(3.21) 



Notice that for the interesting values 5 = 1/2, 1,3/2 we can simplify the resulting expres- 
sions and find 



2' 

5 = 1, 

3 

S= 2' 



E(N) 



N 



+ 1 



V(l) 



E{N) = 4 [4>{N + 3) - ^(3)] 



E(N) 



N 



N 



iP - + 3 +<M- + 2 -^(3) -^(2) 



(3.22) 
(3.23) 
(3.24) 
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The paired states with q 3 = ±q* are associated with degenerate values of the energy. Some 
of them can have zero momentum. 

An example is the case N = 6 at s = \. The Baxter function associated with q 3 = is 

. R 19u 4 323u 2 153 22 , P 

Q(u)=u 6 + , 7=—, e lP = l. (3.25) 

^ K ' 4 80 320 3 

There are other two solutions with P = which are obtained with q$ = ±2 \/723. The 
associated Baxter functions are related by parity. One of them reads 



, fi 2-V/723W 5 235w 4 5 /24T , 2523u 2 23 V~ U 155 

P(u) = ii 6 + — 1 1 \ u 3 = 1 , (3.26) 

^ ; 13 52 143 V 3 2288 1144 9152' V ' 

and has 

7 = f, *" = !. (3.27) 

odd JV 

For even iV, the polynomial P is even 

^(33) = ^(-93)- (3-28) 

Hence we have only q 3 = ±q* associated with degenerate states with Baxter functions 
related by u — > —u. Again, some of these states can have zero momentum. 

An example is the case N = 3 at s = \. There are two paired zero momentum states 
with q 3 = ±§ -JZS and Baxter function (the other is related by parity) 

Q( U ) = u 3 3 /5 2 u _ 1 15 tP 

^ V ' 2 V 7 4 8A/35 2 ' v ; 



3.3 Twist 4 

For L = 4 and spin i\T the Baxter equation is 

[u + t sf Q{u + i) + {u-i sf Q{u -i) = t(u) Q(u), (3.30) 

with 

t{u) = 2 u i + q 2 u 2 + q z u + qi, (3.31) 

q 2 = N - N 2 - &Ns - 12s 2 . (3.32) 

In this case, there are two quantum numbers. They must agree with the multiplicities in 

[s] ® [s] ® [s] ® [s] = [4s + m + n 2 + n 3 ] = ^ + ^ + 2) [4s + (3-33) 

If N is even, looking for a polynomial solution to the Baxter equation we find the condi- 
tions 

P(q 3 , qi ) = 0, (3.34) 
q 3 R{q 3 , qi ) = 0. (3.35) 
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If q 3 = 0, we find 



P(0, g 4 ) = 5( 94 ) = 0, degS = y + 1- (3-36) 



These are non-degenerate states with Q(u) = Q(—u) hence zero momentum. Notice that 
in this case the transfer matrix is even, a property which is related to the parity invariance 
of Q. 

The solutions with q 3 ^ appear in degenerate pairs and can have zero momen- 
tum. If N is odd, looking for a polynomial solution to the Baxter equation we find again 
conditions 

P(q 3 ,q 4 ) = 0, (3.37) 
qsR{qs,qi) = 0. (3.38) 

If q 3 = 0, we find solutions with Q(u) = —Q(—u) hence non-zero momentum. The solu- 
tions with qz ^ appear in degenerate pairs and can have zero momentum. 

3.4 Twist 5 

For L = 5 and spin N the Baxter equation is 

(u + i sf Q(u + i) + [u-i sf Q(u - i) = t(u) Q{u), (3.39) 

with 

t(u) = 2 u 5 + q 2 u 3 + q z u 2 + <? 4 u + q 5 , (3.40) 

q 2 = N - N 2 - 10 Ns - 20 s 2 . (3.41) 

In this case, there are three quantum numbers. They must agree with the multiplicities in 
[.]- = © (g+1>(W + 2>(W + 3) p. + AT], (3.42) 

N=0 ^ 

We focus on the non-degenerate states. These are present for even N and when the trans- 
fer matrix has definite parity. In this case, this means 

t(u) = 2 u 5 + q 2 u 3 + q A u. (3.43) 

The Baxter equation reduces to a polynomial in g 4 that turns out to have degree 

P(g 4 ) = 0, degP = y + l, (3.44) 

as in the L = 4 case. 
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4. Lessons from the previous analysis and general features 



Let us stop to illustrate a few important features emerging from the previous long and 
explicit discussion. 

1. The case L = 2 is well known. There is a single st(2) highest state for each spin. It 
has the correct zero momentum only for even spin. 

2. The case L = 3 is also rather well known [|||, [15|, [16|, |7[ fL8|]. We consider only 
the zero momentum states which are the only relevant ones to planar M = 4 SYM. 
For even spin, the minimal energy state is a singlet. The other excited states appear 
always in degenerate pairs. For odd spin, there are no singlet states. Degenerate 
states are associated with a symmetry in the planar limit relating traces with re- 



versed traces [30] 



Tr(D ni <p---D nL <p) Tr(D nL <p---D ni <p). (4.1) 

3. For L > 3, the (zero momentum) highest weight states can be divided into two 
subsets. Singlets with non degenerate energy, and Paired states with degeneracy 2. 

From the symmetry of the Baxter equation it is easy to proof the 

Theorem 4.1 The singlets are all obtained by solving the Baxter equation with the requirement 
that the transfer matrix eigenvalue t(u) has definite parity 

t(-u) = (-l) L t(u). (4.2) 

This sets to zero several quantum numbers. Let us relabel the remaining free quantum 
numbers (conserved charges) as z t . The pattern is clear from the following list of definite 
parity transfer matrices (remember that q 2 is known) 

i 3 (tt) = 2u z + q 2 u, (4.3) 

U{u) = 2 u 4 + q 2 u 2 + zi, (4.4) 

t 5 (u) = 2 u 5 + q 2 u 3 + zi u, (4.5) 

t % {u) = 2u 6 + q 2 u i + z 1 u 2 + z 2 , (4.6) 

t 7 (u) = 2u 7 + q 2 u 5 + ziu 3 + z 2 u. (4.7) 

The number of singlet states for a certain twist is the number of possible values of these 
quantum numbers. It is a function of the Lorentz spin given by the following simple 
formula 

L = 2n,2n+1, #singlets= 2 ^ . (4.8) 

For illustration, we show in Figs. ([I]) and (||) the full spectrum of highest weight states 
at L = 3, 4. A general feature is that the singlet part of the spectrum embraces the full 
spectrum. In particular, the lowest and highest states are singlets. 
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5. Logarithmic scaling of the anomalous dimensions 



The following general information is known about the band of highest weight anomalous 
dimensions at generic twist L. We focus on the scalar s = 1/2 sector, but generalizations 
are possible. In the N — > oo limit, the minimal anomalous dimension has a logarithmic 
scaling which is twist independent and reads 

7min ~ f(g) lOg N. (5.1) 

The most explicit proofs of this statement in M = 4 SYM are [[31, |2|, |33|] and [35| ]. 
The scaling function f(g) is proportional to the physical coupling, a.k.a. cusp anomalous 
dimension [ j3~ll ] and has been computed at all orders in p^ ]. At three loops, it reads 

/(3)=4 5 ^ h , (5.2) 

2 7T 2 4 11 7T 4 fi ._ 

= 9 f o H . (5.3) 

y 6 y 180 a 



In twist-2, the physical principle behind the scaling Eq. fl5.1|) is simply that the large N 
limit is nothing but the quasi-elastic aiBjorken — > 1 deep inelastic scattering regime. This 
is dominated by universal classical soft gluon emission characterized by the anomalous 
cusp contribution of the quasi-free parton dynamics. 

The excited anomalous dimensions are expected to scale in the same way but with a 
different prefactor ranging as follows 

f(g) log N < 7 < \ f(g) log N. (5.4) 

We remark that this can be nicely understood, at strong coupling, in terms of the dual 
string configurations which have L spikes each contributing | f(g) to the coefficient in 
the case where they are equally spaced [|24|]. 

At next-to-leading logarithmic order, a general formula has been recently derived 



in pq ] for the minimal anomalous dimension at generic twist. It reads 

7min = f(g) log N + f s i(g, L) + suppressed terms, (5.5) 
where, in our notation for the coupling 

M(9, L) = (7s - {L - 2) log 2) f(g) - 2 (7 - 2 L) C 3 <? 4 + (5-6) 
+ 7T 2 C 3 + (62- 21 L)( 5 ) 5 6 + -- - . 



This formula is remarkable because it gives the explicit twist dependence, reabsorb the 
scaling function in a compact way and provides the other corrections as £-terms with 
simple twist dependence. 

Given the above general constraints, is it possible to explore analytically the spectrum 
of highest weights at general twist ? A partially positive answer is provided by sum rules 
that we know describe. 
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6. Linear sum rules at one-loop 



Let 7i j.(iV) denote the anomalous dimensions of the various highest states, labeled by k. 
We consider at one-loop all the three physical values s = 1/2,1,3/2 associated with ele- 
mentary scalars, gauginos, and gauge fields. We can compute the sum of the anomalous 
dimensions of singlet states 

4 S V)= E -rSW- (6-1) 

fc£ singlets 

It turns out that this quantity is rational ! The reason is very simple. The above sum can 
be computed by the Baxter approach. The anomalous dimensions are given by a rational 
function of the free charges which are not fixed by the parity constraint on the transfer 
matrix. These charges are constrained and fully determined by a system of polynomial 
equations. So, the point is to show that the sum of a rational function over the roots of a 
system of (rational) polynomials is rational. This simple theorem is proved and discussed 
in App. (§). 

Given a sequence of rational numbers describing the N dependence of H L (N), it is 
possible to look for closed formulas, by using some trial and error combination of har- 
monic sums, as inspired by the low L cases. Our notation for harmonic sums is standard 
(o6Z,a = (oi,..., a n ) with £ Z) 

U^jffif. S^N ) = Z^fs blk) . (6.2) 

The following closed formulae are obtained for the cases L = 4, 5, 6, 7 (they are fulfilled 
by any N we have been able to test, typically of the order 0(100)) 

K 

L = 4, 5 s£ V) = E [4, 5) in) - a[ 5) (0)] , (6.3) 

n=l 
jv 

L = 6, 7 E«(JV) = E E V) " 4, m) (0)] , (6-4) 

n—l m—1 

where for L = 4 

6Si(n) + 2Si(2n- 1) -2 5_i(2n- 1), (6.5) 
2 5i(n) + 45i(n+ 1) + 25i(2n+ 1) -2S i(2n+ 1), (6.6) 
2 5i(n + 1) +45i(n + 2) + 25i(2n+ 1) -25_i(2n + 1), (6.7) 

for 1 = 5 

8 5i(n), (6.8) 

6 51(71 + 1) + 2Si(2n + 1) - 2 5_i(2n+ 1), (6.9) 

2 5i(n + 1) + 6 5i(n + 2), (6.10) 



°4 



(n) 



r (3/2) 



(n) 



r (V2) 



(n) 

4 3/2) (n) 
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for L = 6 



<4 1/2) (n) = 10 5i(n) + 25i(2n- 1) -25_i(2n- 1), (6.11) 

a^\n) = 2Si(n) + 8Si(n+l) + 2Si(2n+l)-2S_i(2n + l), (6.12) 
4 3/2) (n) = 25i(n + l) + 85i(n + 2) + 25i(2n+l)-25_i(2n + l), (6.13) 

and for L = 7 

4 1/2) (n) = 12 5i(n), (6.14) 

4 1} (n) = 10 5i(n + l) + 25i(2n+l)-25_i(2n + l), (6.15) 

4 3/2) (n) = 2 5i(n+l) + 105i(n + 2). (6.16) 

Notice that for s = 1/2 we can also write in a more uniform way 

4 1/2) (n) = 8 5i(2n) + 4 5_i(2n), (6.17) 

4 1/2) (n) = SSxin), (6.18) 

4 1/2) (n) = 12 Si (2 n) + 8 5_i(2n), (6.19) 

4 V2) (n) = 12 5i(n). (6.20) 

where we have exploited the remarkable identity 

5i(2s-l)-5_i(2s-l) = 25_i(2s)+45i(2a)-35i(s), s g 2N. (6.21) 



7. Linear sum rules at three loop results in the scalar sector 

Starting from the one-loop Bethe roots evaluated with the Baxter approach, one can build 
the multi-loop anomalous dimensions by feeding the long-range Bethe equations of 
This is quite easy in the s = 1/2 case where the Bethe equations are not nested. The 
procedure is standard (see for instance the detailed discussion in [[l5|]). From a long list 
for several even values of N, one makes an Ansatz with higher transcendentality nested 
sums and solves the over constrained system of equations. Dropping for simplicity the 
s = l/2 label and denoting 

^(n) = ^ 5 2 V v (n), (7.1) 

one finds the following solutions. 
L=4 

The argument of the harmonic sums is 



S... = S...{2n). 



(7.2) 



<T4,2 = 16 5_3 + 24 53 — 16 S_2,-i — 8 5_2,i — 8 5_i i _2 — 8 5_i, 2 — 16 5i,_ 2 

-16 51,2 - 16 52-1 - 24 52,1, (7.3) 
^ 3 /2) = 104 5_ 5 + 184 5 5 - 144 5_ 4 ,-i - 104 5_ 4 ,i - 192 5_ 3 ,- 2 - H2 5_ 3 , 2 - 144 5_ 2 ,- 3 

-96 5_ 2 ,3 - 96 5_i _ 4 - 96 5_i, 4 - 192 5 X _ 4 - 192 5i, 4 - 184 5 2 ,_ 3 - 248 5 2 , 3 - 176 5 3 _ 2 
-288 5 3 , 2 - 144 5 4 ,_i - 216 5 4 ,i + 64 5_ 3 ,-i,_i + 128 5_ 3 ,-i,i + 64 5_ 3 ,i,_i + 32 5_ 3 ,i,i 
+64 5_ 2 ,_ 2 ,_i + 96 5_ 2 ,_2,i + 64 5_ 2 ,-i,-2 + 64 5_ 2 ,-i,2 + 32 5_ 2 ,i,-2 + 32 5_ 2 ,i,2 
+64 5_ 2 ,2,_i + 32 5_ 2 ,2,i + 64 5_i,_ 3 ,_i + 64 5_i,_ 3 ,i + 32 5_i,_ 2 ,_ 2 + 32 5_i,_ 2 , 2 
+32 5_i, 2 ,-2 + 32 5_i, 2 ,2 + 64 S-i, 3 ,-i + 64 5_i, 3 ,i + 128 Si,- 3 ,-i + 128 Si,_ 3 ,i + 64 5i,_ 2 ,-2 
+64 5i,_ 2 ,2 + 64 5i, 2 ,-2 + 64 5i, 2 , 2 + 128 5i, 3 ,_i + 128 5i, 3 ,i + 128 5 2 ,_ 2 ,_i + 80 5 2 ,- 2 ,i 
+64 5 2 ,_i,_ 2 + 64 5 2 ,_i, 2 + 96 5 2 ,i,_ 2 + 96 5 2 ,i, 2 + 128 5 2 , 2 ,-i + 160 5 2 , 2 ,i + 128 S 3 ,-i,-i 
+64 5 3 ,_i,i + 128 5 3 ,i,_i + 192 5 3 ,i,i . (7.4) 



L=5 

The argument of the harmonic sums is 

5... = 5...(n). (7.5) 

£7 5 ,2 = 8 5 3 - 8 5i, 2 - 12 5 2 ,i, (7.6) 
a 5 , 4 = 21 5 5 - 24 5i, 4 - 38 5 2 , 3 - 46 5 3 , 2 - 36 5 4 ,i + 16 5i, 2 , 2 + 32 5i, 3 ,i 

+24 5 2 ,i, 2 + 40 5 2 , 2 ,i + 48 5 3 ,i,i . (7.7) 

L=6 

The argument of the harmonic sums is 

5... = 5...(2n). (7.8) 

C7 6 , 2 = 40 5_ 3 + 48 5 3 - 32 5_ 2 ,_i - 24 5_ 2 ,i - 16 5_i,_ 2 - 16 5_i, 2 - 24 5i,_ 2 

-24 5i, 2 - 32 5 2 ,_i -40 5 2 ,i, (7.9) 

ct 6i 3 = 576 5_ 5 + 736 5 5 - 480 5_ 4 ,_i - 384 5_ 4 ,i - 560 5_ 3 _ 2 - 448 5_ 3 , 2 - 400 5_ 2 ,-3 
-352 5_ 2 ,3 - 192 S_i _ 4 - 192 5_i, 4 - 288 5 X _ 4 - 288 5i, 4 - 464 5 2 ,_ 3 - 512 5 2 , 3 
-560 5 3 ,- 2 - 672 5 3 , 2 - 480 5 4 ,_i - 576 5 4 ,i + 256 5_ 3 ,_i,_i + 320 5_ 3 ,_i,i 
+256 5_ 3 ,i,-i + 192 5_ 3 ,i,i + 192 5_ 2 ,-2,-i + 224 5_ 2 ,_ 2il + 128 5_ 2 ,-i,-2 + 128 5_ 2 ,-i,2 
+96 5_ 2 ,i,-2 + 96 5_ 2 ,i,2 + 192 5_ 2 , 2 ,_i + 160 5_ 2 ,2,i + 128 5_i,_ 3 ,_i + 128 5_i,_ 3 ,i 
+64 5_i,_ 2 ,-2 + 64 5_i,_ 2 ,2 + 64 5_i, 2 ,_ 2 + 64 5_i, 2 , 2 + 128 S_i, 3 ,-i + 128 S_i, 3 ,i 
+192 5i,_ 3 ,_i + 192 5i,_ 3 ,i + 96 5i,_ 2 ,_ 2 + 96 5i,_ 2 , 2 + 96 5i, 2 ,- 2 + 96 5i, 2 , 2 + 192 5i, 3 ,_i 
+192 5i, 3 ,i + 256 5 2 ,_ 2 ,_ 1 + 224 5 2 ,_ 2 ,i + 128 5 2 ,_i,_ 2 + 128 5 2 ,-i, 2 + 160 5 2 ,i,_ 2 
+160 5 2 ,i, 2 + 256 5 2 ,2,-i + 288 5 2 , 2 ,i + 320 5 3 ,_i,_i + 256 5 3 ,_i,i 

+320 53,1,-1 + 384 53,1,1. (7.10) 
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L=7 

The argument of the harmonic sums is 

5... = 



S...(n). 



(7.11) 



a 7;2 = 14 53-12 51,2-20 52,1, 

(7 7 ,3 = 61 5 5 - 36 5i, 4 - 70 5 2 , 3 - 94 5 3 , 2 - 84 5 4 ,i + 24 5i, 2 , 2 + 48 5i, 3 ,i 
+40 5 2 ,i,2 + 72 52,2,i + 96 5 3 ,i,i . 



(7.12) 



(7.13) 



Notice that for s = 1/2 the general expressions 



L = 



4,5 



E[4 s) N-4 s) (o)] 



(7.14) 



n=l 



N 



L = 6,7 



4 s) w = E E [4 s V)-4 m) (o)] 



(7.15) 



n=l m=l 



simplify since a^ 2 ^ (0) for all the considered L and up to 3 loops. 



8. Linear sum rules: Structural properties and twist-dependent formulas 

The previous results show various remarkable structural properties. These are 
1. The general formula for Sj, up to three loops takes the form 



where the number of sums is p = n — 1 for both L = 2 n and L = 2 n + 1. 

2. The internal function a^rip) can be written as a linear combination of harmonic 
sums with total transcendentality equal to 2 1 — 1 where I is the loop order i = 1,2,3. 

3. The argument of the harmonic sums is n p for odd L and 2 n p for even L. 

4. The multi-index of the harmonic sums does involve only positive indices for odd L. 

5. The set of multi-indices is the same for all even L and fixed loop order. The same is 
true for odd L with a different set of indices. 

We have extended the calculation up to L = 13 testing the above structural properties. In 
all the considered cases they hold. Also, looking at the L dependence of the coefficients of 
the harmonic sums, we have been able to write down the following compact and, in our 
opinion, remarkable expressions. 



El(W) 



E E ■ ■ ■ E a L{np) 



(8-1) 
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8.1 Odd twist 

Up to three loops, we have 

X L {N)= 2(L-l)S x , ig 2 + (8.2) 
+ [(3 L - 7) S x , 3 — 2 (L — 1) 5 x ,i,2 — 4 (L — 2) 5 x , 2 ,i] 3 4 + 
+ [(20 L - 79) 5 X , 5 — 6 (L — 1) S XAA - 12 (2 L - 7) 5 x , 4 ,i + 
-2 (8 L - 21) 5 X ,2,3 — 2 (12 L — 37) 5 X , 3 ,2 + 4 (L — 1) 5 x ,i,2,2 + 
+8 (L - 2) Sjc i2i i i2 + 8 (2 L - 5) S x , 2 ,2,i + 8 (L - 1) 5 Xil , 3 ,i + 
+24(L-3)5 x ,3,i,i]3 6 + --- • 

where 

5x,a = S x>a (jj, X = {CV-^0}, (8.3) 

L-3 
2 

and a harmonic sum with trailing indices is defined as 

N 

S 0l a(tf) = E 5 ^(n). (8.4) 

This is the usual definition provided we define sign(0) = 1. 

The above formula works in all the considered case. It also works for L = 3 up to two 
loops. The three loop term is not covered for this initial value. 

8.2 Even twist 

Due to the larger computational complexity of the even twist case, we only present a result 
at the two loop level. We define in this case 

5x,a = Sx,a(fV X = {0,---,0} (8.5) 



and (notice the most inner 2 i p argument) 

So,...,o, a (n) = E E £ ■ ■ ■ E S a (2i p ) (8.6) 



One finds for even L > 4 



X { l /2 \N) =[2L S x ,i + 2 (L - 2) S x ,-i] Q 2 + (8.7) 
+ [4(3L - 8) S x ,-3 + 12(L - 2) S x>3 - 8(L - 2) S x ,-2,-i ~ 8(L - 3) S*,_ 2 ,i 
-4(L - 2) 5 x ,-i,-2 - 4(L - 2) 5 x ,-i, 2 - 4L 5 x ,i,- 2 - 4L 5 x ,i, 2 
-8(L - 2) 5 Xi2i _i - 8(L - 1) 5 x ,2,i] 3 4 + • • • • 
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9. Large N check: Recovering the cusp anomalous dimension 

An important check of the previous results Eqs. ( |&2| , |8^ ) is that for large N all the ground 
and excited anomalous dimensions are expected to scale logarithmically with N with a 
coupling dependence reabsorbed in the physical coupling g^ h . 

To check this, let a = (ai , a 2 , ■ ■ ■ ) and = . If ai 7^ 1 we have at large N 

Sx,id N ) ~ — 5 a(oo) log N, X = 0^_0 . (9.1) 
V ' v 
Apart from trivial factors depending on the multiplicity of the singlet set of states, we can 
read the coefficient of the logarithmic leading term by collecting all the nested harmonic 
sums with a leading 1 index and replacing 

Si,s{N) -> S s (oo). (9.2) 

We now present the detailed check of the g^ h reshuffling for all the expressions that we 
have listed in the previous sections. 

9.1 Odd twist at three loops 

In this case, we use Eq. d&2| ) and write 



El(#) ~ A ( ?V K9) log N, P = ^, (9.3) 



p! 

where the function h(g) is 

h(g) = 2{L-l)g 2 + (9.4) 
+ [-2(L-l)5 2 (oo)] 5 4 + 

+ [ - 6 (L - 1) S 4 (oo) + 4 (L - 1) 5 2 , 2 (oo) + 8 (L - 1) 5 3 ,i(oo)] 5 6 + • • • . 
Notice that the number of singlets is asymptotically 

# S ingle, S =(«;*)~i(f) r (9.5) 

Hence, we can divide by the multiplet dimension and write 

E L (i\T) ~ h{g) log N. (9.6) 
Replacing the following asymptotic sums 



S 2 (oo) = C2 = y> 5 2, 2 (oo) = — , 



(9.7) 



4 4 

7T „ , , 7T 



we find 



5 4 (oo) = £ 4 = — , 5 3 ,i(oo) 72 



%) = 2(L-1) ^_ C25 4 + 1^! 5 6 + 



2(L-l) S ; h . (9.8) 
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9.2 L = 4 at three loops 

The quantities cr 4 ^ contain the following nested sums with a leading 1 (and argument 2 n) 



(74,1 = 8Si + --- , (9.9) 

<7 4 ,2 = -16(5i, 2 + 5i_ 2 ) + --- , (9.10) 
(7 4 ,3 = -192 (Si, 4 + 5i _ 4 ) + 128 (5i, 3 ,i + 5i, 3l -i + Si,- 3 ,i + Si,- 3 ,-i) + 

+64 (Si.2,2 + Si.2,-2 + 5i_ 2l 2 + Si-2,-2) + " " " ■ (9-11) 

The required asymptotic values are 

5 2 + S_ 2 | 00 = ^, (9.12) 

s « + s -«loo = ^. ( 9 - 13 ) 

5 3 ,i + S 3l _i + S_ 3 ,i + ^3,-iL = ^g, (9.14) 

S 2> 2 + S 2 ,_ 2 + S_ 2 , 2 + = (9.15) 



Collecting, we find 



a 4 ~ 8 log N (V -C25 4 + ^-5 6 + ---j = 

= 8 1ogiV 5 2 h . (9.16) 

9.3 L = 6 at three loops 

The quantities contain the following nested sums with a leading 1 (and argument 2 n) 

(7 6 ,i = 12 5i + --- , (9.17) 

(7 6 , 2 = -24 (Si i2 + Si,_ 2 ) + --- , (9.18) 
a 6 , 3 = -288 (5i, 4 + 5i_ 4 ) + 192 (5i, 3 ,i + Si, 3 ,_i + Si_ 3 ,i + Si-3,-1) + 

+96 (5i i2 ,2 + Si i2 ,_2 + 5i,_ 2i2 + Si,_ 2i _ 2 ) + • • • . (9.19) 

With the previous asymptotic values, we find 

C7 6 ~ 12 log N - ( 2 5 4 + g 6 + • • • j = 

= 12 1og (9.20) 

9.4 General even L at two loops 

The asymptotic values given for the cases L = 4, 6 are sufficient to check the two loop 
general even L case. 
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10. Linear sum rules: Subleading corrections at large N 



We can take our master formula for odd twist Eq. ( |&2| ) and expand Sl(2V) at large N com- 
puting the subleading part analogous to / s i defined in Eq. The expansion requires a 
treatment of multiple sums with various trailing zeroes of the form 



So,...,o,x- 



(10.1) 



They can be treated with the nice results reported in App. Q). It is clear that some struc- 
ture immediately arises. For instance, let us consider the twist 5 case where p = 1. The 
leading terms in S , proportional to N, arise from sums of the form So,i,a,x(rc) where 
n = N/2 and a > 1. From Eq. ( |C.5[ ) of App. Q, we get 



Sb,i,a,x(n) = (n + 1) Si, a ,x(n) - 5 , a ,x(n) = 

= (n + 1) [5i i0 , x (n) - 5 , x (n)] + S a -i,x(n). 

For large n, we need only the first bracket whose expansion contains the terms 

Si,a,x(ra) - 5 a ,x(n) = (log n + 7 B - 1) £a,x + . . . , 



(10.2) 
(10.3) 



(10.4) 



where dots denote other constant terms. The first part combines to give the physical 
coupling which thus must appear in the combination 



(log n + 7s - l)gp h . 



(10.5) 



The prefactor is easily generalized to a generic odd twist, i.e. p > 1. The leading term of 
So,...,o,i,x{n) for large n replaces it by the general form 



log n + 7e - 1 



p\ 



[a ?; (logn + 7E) - b p ] 



(10.6) 



where 



dp — O/p— i — lj 

bp = , bi = 1, — > I 

P 

The remaining subleading pieces can be worked out in a similar way. The final result is 



(10.7) 
(10.8) 



E L (tf) = 2(L-1) 



',2V „ /L-3 

log y +7E-M^ 



3ph + 



(10.9) 



(3L-7)Cs5 4 + 



+ 



L-2 



7T 2 Cs + (lOi - 31) Cb 



5° + 



Eq. fll0.9j ) is the generalization of the recent result Eq. (^). The structure is quite similar, 
although Eq. ( |10.9j ) involves the sum over the singlet anomalous dimensions ! 
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11. Higher order sum rules: The quadratic case 

Let us define higher order sum rules by considering sums of powers of the individual 
anomalous dimensions. In particular, we focus on the quadratic sum 

Q { L ) m= E bKw] 2 - (ii.i) 

fc£ singlets 

Remarkably, we find simple sum rules also for these higher order sums. We illustrate this 
in the special cases L = 4, 5 for s = 1/2. We find again the general representation (valid 
for L = 4,5) 

QL(N) = J2E9 2i+2 ^An), (H.2) 

n=l 1>1 

where 
L=4 

The argument of the harmonic sums is 

S... = S...(2n). (11.3) 



g 4 ,i = -48 5_2 - 128 5 2 + 64 5_i_i + 32 S_ M + 64 5i _i + 128 Si,!, (H-4) 
94,2 = -384 S_ 4 - 832 S 4 + 576 S_ 3 ,_i + 256 S_ 3 ,i + 704 S_ 2 _ 2 + 384 5_ 2 , 2 

+576 S_i ,_3 + 384 S_i,3 + 768 Si ,_ 3 + 1152 Si, 3 + 768 5 2 ,_ 2 + 1280 5 2 , 2 + 576 5 3 ,_i 

+ 1024 5 3 ,i - 256 5_ 2 ,_i,_i - 512 S_ 2 ,_i,i - 256 5_ 2 ,i,_i - 128 S_ 2 ,i,i - 256 5_i,_ 2 ,_i 

-384 S_i,_ 2 ,i - 256 5_i,_i,_ 2 - 256 5_i,_i, 2 - 128 S_i,i,_ 2 - 128 S_i,i, 2 - 256 S_i, 2 ,_i 

-128 S_i i2i i - 512 Si,_ 2 ,_i - 256 5i,_ 2 ,i - 256 Si_i_ 2 - 256 5i,_i, 2 

-512 Si,i,_ 2 - 512 Si,i, 2 - 512 Si, 2 ,-i - 768 Si, 2 ,i - 512 S 2 ,-i,-i - 256 5 2 ,_i,i 

-512 5 2 ,i,-i - 896S 2 ,i,i. (11.5) 



The argument of the harmonic sums is 

5... = S...(n). (11.6) 



95,1 
95,2 

95,3 



128 5i,i - 88 5 2 , (11.7) 
-212 5 4 + 384 5i, 3 + 432 5 2 , 2 + 352 5 3 ,i - 256 5i,i, 2 

-384 51,2,1 - 448 52,1,1, (11.8) 

-702 5 6 + 1616 Si, 5 + 2112 5 2 , 4 + 2352 5 3 , 3 + 2072 5 4 , 2 + 1408 5 5 ,i 

-1536 5i,i,4 - 2080 5i,2, 3 - 2464 5i, 3 , 2 - 2304 5i, 4 ,i - 2336 5 2 ,i, 3 

-2784 5 2 , 2 , 2 - 2768 5 2 , 3 ,i - 2912 5 3 ,i, 2 - 2912 S 3 , 2 ,i - 2544 5 4 ,i,i + 768 5i,i, 2 , 2 

+ 1536 5i,i, 3 ,i + H52 5i, 2 ,i,2 + 1920 Si, 2 , 2 ,i + 2304 Si, 3 ,i,i + 1344 5 2 ,i,i, 2 

+2112 5 2 ,i,2,i + 2496 5 2 ,2,i,i + 2688 S 3 ,i,i,i. (11.9) 
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12. Quadratic sum rules: Structural properties and twist dependent formulas 

In complete analogy with the linear sum rules, we can observe basically the same struc- 
tural properties also in the case of the quadratic sum rules. In particular the general for- 
mula for Q L seems to be 

Ql(N)= f: E ■■■ E QlM (12.1) 

ni = l «2 = 1 n v =l 

where the number of sums is again p = n — 1 for both L = 2 n and L = 2 n+ 1, but now the 
total transcendentality of the sums in is equal to 2 £ where £ is the loop order £ = 1,2,3. 

Also in this case, we have extended the calculation up to L = 13. Now, the L depen- 
dence of the harmonic sums coefficients is quadratic instead of linear and we can write 
down the following compact expressions. 

12.1 Odd twist 

We have, at two loops, 

Q L (N) = [8 (L - l) 2 S*,i,i - 4 (2L 2 -1L + l) S x<2 \ g A + 

+ [ — 16 (L — l) 2 Sjr.1,1,2 + 12 (3L - 7)(L - 1) S XX3 
-32 (L - 2)(L - 1) 5 x ,i,2,i - 4 (l2L 2 - 74L + 123) 
+4 (l3L 2 - 58L + 73) S x ,2,2 + 8 (7L 2 - 38L + 59) 
-16 (3L 2 - 12L + 13) S x ,2,i,i] 9 6 + • • • • 

where, again, 

Sx,, = S x , a (y) , X = {^^0} 

L-3 
2 

The cusp anomaly check is clearly passed by the combination 
8 (L - l) 2 Sjc.i.i 5 4 — 16 (L — l) 2 Sjc.1,1,2 5 6 . 

12.2 Even twist 

Due to the larger computational complexity of the even twist case, we only present a one- 
loop result. We define in this case 

Sx, a = S x>a X = {CV^O} (12.5) 

L-2 
2 

and (notice the most inner 2 i p argument) 



(12.2) 



Sx,i 

Sx,3,l 



(12.3) 



(12.4) 



so o, a (n) = J2 EE-EW 



(12.6) 



One finds for even L > 4 

Ql(N) = [ - 8 (2L 2 - 5L + 4) S xa - 8 (2L - 5) (L - 2) 5 X ,_ 2 
+8L 2 5 X , 1 , 1 + 8L(L-2)5 X , 1 ,_ 1 

+8 (L - 2) 2 5 x ,-i,i + 8 L (L - 2) Sjr._i._i] <? 4 + • • • (12.7) 
13. Large N check for the quadratic sum rules 

Here we report the cusp anomalous dimension check for the various formulas computing 
quadratic sum rules at L = 4 (two loops) and L = 5 (three loops). 

13.1 L=4 at two loops 

The squared logarithmic terms are 

g 4l i = 1285i,i + ■ ■ ■ , (13.1) 

94,2 = -512Si.i _ 2 - 5125i,i.2 + • • • . (13.2) 

Hence, at two loops 

94 = 64 log 2 N(g* -2C2 9 6 + ■■■) = 

= 64 log 2 iV( 5 2 -C 2 5 4 + ---) 2 = 

= 64 1og 2 iV«4. (13.3) 

13.2 L=5 at three loops 

The squared logarithmic terms are 

95.1 = 1285i.i + ■ ■ ■ , (13.4) 

95.2 = -25651,1.2 + • • • , (13.5) 

95.3 = -15365i,i, 4 + 7685i,i i2 ,2 + 15365i,i, 3 ,i + ■ ■ ■ . (13.6) 

Hence, at three loops 

95 = 64 log 2 N ( 5 4 - 2 c 2 g 6 + ^ g 8 + ■ ■ ■ ) = 

■ -.4 1 11 „6 1 ^2 



64 


log 2 


N(g* 


64 


log 2 


N(g 2 


64 


log 2 





14. One loop cubic sum rule 



(13.7) 



It is clear that it is possible to derive sum rules at arbitrary high order. The determination 
of the explicit formulas is a matter of computational effort. Here, we just give, as an 
example, the cubic sum rule 

fc£ singlets 
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for the scalar sector s = 1/2 and odd twist L at one-loop. We have tested it again up to 
L = 13. It reads 

C L (N) = [48 (L - l) 3 S XAAA — 24 (L — 1) (2L 2 -7L + T) S x>1>2 
-48 (L - 2) (l, 2 - 4L + 7) 5 x , 2il 

+8 (6L 3 - 45L 2 + 124L - 12l) S Xfi ] g 6 + ■ ■ ■ . (14.2) 



where, again, 



Sx, a = %,a ( f ) , X = {0,--- ( 0} (14.3) 



Z/-3 
2 



15. Conclusions 

In summary, we have shown that it is useful to define higher-order sum rules for the 
anomalous dimensions of singlet unpaired twist-operators in JV = 4 SYM for arbitrar- 
ily high twist. Of course, these combinations contain less information than the separate 
anomalous dimensions. However, on the other hand, they admit multi-loop closed ex- 
pressions in terms of the usual nested harmonic sums which are ubiquitous in this con- 
text. These expressions provide interesting handles for analytical calculations in higher 
twist. 

The basic technical hint behind the sum rules is that they are related to restricted 
traces of powers of the dilatation operator. As such, they are quite simpler objects than 
the separate anomalous dimensions. As an analogy, it is typically simpler to focus on the 
coefficient of a polynomial instead of looking at its explicit and complicated roots. 

Our multi-loop results have been obtained in the sl(2) sector where the long-range 
Bethe equations are particularly simple. At one-loop, analogous results have been pre- 
sented for the other basic sectors describing purely fermionic or gauge operators. It 
should be possible to extend the analysis to these cases at higher loops by working out 
the perturbative expansion of the relevant higher rank long-range equations. 

It remains to be understood if the closed expressions we found for the sum rules are 
just a curiosity or a manifestation of deeper properties. In this respect, their interpretation 
in the light of AdS/CFT duality would certainly be a very interesting issue. Indeed, in this 
context, the linear sum rules compute sums of energies of dual string configurations with 
a fixed number of spikes, but different values of the internal degrees of freedom associated 



with the band of states (for twist > 2) [p2|]. Hence, on the string side, the proposed sum 
rule suggest to investigate the properties of energies after a sort of averaging over these 
kinematical features. The powerful analysis in [24 1 could prove to be useful in this respect. 
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A. Symmetric polynomials and sums of powers 

Given a finite set of complex numbers x\ , . . . , x n £ C, the symmetric polynomials life [x\ , . . . , x n ) 
are defined as 

f\(x - Xl ) = n fc (xi, . . . , X n ) X k . (A.l) 

i=l k=0 



Hence, 



U {x 1 ,...,x n ) = 1, (A.2) 

Ili(a!i, . . . , x n ) = X\ H x n , (A.3) 

n 2 (a;i, ■ ■ ■ ,x n ) = XiX i> ( A - 4 ) 



l<i<j<n 



(A.5) 

n 

n n (an, . . . ,x n ) = Yl x i- ( A - 6 ) 



1=1 

The relation with the symmetric sums of powers 



k 



i=l 

is given by the generating function relation 



^n fc i fc = exp £(-l) fc+1 ^i fc , (A.8) 
k=0 \fc=i K J 

or by the Newton-Girard recursion 

m 

(-irmn m + ^(-l) fc+m n m _ fc 5 fc = 0. (A.9) 
k=l 

The resulting map (111, . . . , n„) <-> (Si, . . . , S„) is birational, actually polynomial, and 
begins with 

Si = Hi, Hi = Si, 

s 2 = n? - 2 n 2) n 2 = \ {s\ - s 2 ) , (A.io) 

s 3 = nf -3n!n 2 + 3n 3l n 3 = £ (s? - 3Si s 2 + 2S 3 ) . 

B. Rationality proof 

Theorem B.l Let R(x) be a rational function of x over Q, i.e. the ratio of two polynomials in 
Q[x]. For any polynomial with rational coefficients P{x) £ Q[x] one has 

J2 R( x ) e Q • (B.l) 

xGC: P(x)=0 
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Proof: 

Let d = deg P. For any root of P(x) = we can write the identity 



d-l 



R{x) = Y J CnX n , (B.2) 

with suitable coefficients {c n } which are rational functions of the coefficients appearing 
in R. They are the same for all roots of P. As is well known (see App. (|A|)), the sums 

S n = E ( R3 ) 

sfEC : P(x)=0 

are all fully determined as rational functions of the coefficients of P(x). Thus they are 
rational and therefore 

d-l 

E = E c„ S n 6 Q. (B.4) 

xGC: P(x)=0 n=0 



This theorem can be greatly extended by replacing a; by a finite set of variables and 
P(x) = by a system of polynomial equations with rational coefficients. The proof is the 
same as in the one dimensional case after reduction by Grobner basis methods [^]. 

In order to illustrate the above (constructive) theorem in the univariate case, let us 
consider in details as a simple example the derivation of 

E i = -l- (B-5) 

x 6 +ir+l=0 

We start from the relation 

- = -l-x 5 , iff z 6 + a: + l = 0. (B.6) 
x 

Hence 

E | = -8- E * 5 - (B-7) 

E 6 +E+1 = X 6 +E+1 = 

From the results of App. @ we compute for P(x) = x 6 + x + 1 the explicit sums of powers 

S 1 = S 2 = S 3 = S i = 0, 5 5 = -5. (B.8) 
Hence, we have proved that 

E i = -6-S B = -l. (B.9) 

E 8 +E+1 = 
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C. Harmonic sums with trailing indices and positive indices 

A leading index means 

N 

So,x(tf ) = E 5 x(n). 



n=l 



Clearly, we have 



S (N) = N, 

N(N + 1) 



and the general formula 



S ,o{N) 
So o(N) 



2! ' 

'N + p- i x 

p 



(C.1) 

(C.2) 
(C.3) 

(C.4) 



The first non trivial result is 
Theorem C.l For any a > 1, we have 



(C.5) 



Proof: 

We can prove the theorem by induction on N. Taking the difference of the equation 
between N + 1 and N we find 

So l0 ,x(tf +l)-(W+2) 5a,x(^+l)+5 a _i,x(^+l)-5o, a ,x(iV)+(iV+l) S a ,x(N)- S a . h x(N) = 
= S a ,*(N + 1) - (N + 2) 5 a , x (iV + 1) + ^-L-^ iV + 1) + (N + 1) S , x (tf ) = 



SaM N + 1) " Sa,x(W) 



= 5 a , x (iV + 1) - 5„,x(JV) 

The generalization is 

Theorem C.2 For any a > 1, we have 



1 



(JV+ 1)' 



S X (W + 1) = 0. 







V^a.xW = - 
p 



(JV + p) 5 0) ... l o,a,X - <So,..,0,a-l,X 



p-l 



p-l 



(C.6) 



Proof: it follows from induction over p. Let us assume that Eq. ( |C.6D is true for p for all N, 
then we can prove that it is true for p + 1 by induction over N. Following similar steps as 
in the proof of Theorem C.l we find 



S ,..., ,a,x( N + 1 ) 



P+l 



p+l 



(N + p + 2)S 0> ... Aa>x (N + 1) - S ,... Aa - lt x(N + 1) 
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So 0,a,x( N ) + 



P+1 



p+1 



(N + p+ l)5 ,..., ,a,x(^) - S , a -l,xW 



So,...,o,a^{N + 1) - N +* + 1 S , a , x {N + 1) - -^S 0t ... Aa>x (N + 1) 



p + 



p-1 



p + 1 



(N + p + l)S ,..., ,a,x(W + 1) - 5 ,...,0,a-l,x(iV + 1) 



p-i 



p-i 
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Figure 1: Full spectrum at twist L = 3. 
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Figure 2: Full spectrum at twist L = 4. 
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